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We give comments on the recent paper by Giunti (Phys. Lett. B 665, 92 (2008), arXiv: 0805.0431 [hep-ph] 
) with a critique of our explanation of the experimentally observed periodic time-dependence of the interference 
term in the rate of the K-shell electron capture decay of the H-like ions 140 p r 58 + and 142 Pm 60+ as a two-neutrino- 
flavour mixing. We show also that this phenomenon cannot be explained by a coherent mixing of two states of a 
mother ion as proposed by Giunti. 
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Recently Litvinov et al. PQ have observed that 
the K-shell electron capture {EC) decay rates of 



the H-like ions 140 Ce 58+ or 142 Nd' 



60+ 



140p r 58+ 
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have unexpected oscillatory structure. According 
to the experimental data [1] , the rates of the num- 
ber Nf c of daughter ions 140 Ce 58+ or 142 Nd 60+ 



dt 



(2) 



where N m (t) is the number of mother H-like ions 

140p r 58+ Qr 142p m 60+ p and A (H) ^ ig ^ E(J _ 

decay rate, are periodic functions, caused by a 
periodic time-dependence of the EC-decay rates 



^(*) = A^(l + 

with a period Td 
0.20. 
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We have proposed in [5] an explanation of the 
periodic time-dependence of the .EC-decay rates 
as the interference of two massive neutrinos v\ 
and i>2 with masses mi and mi , respectively. The 
period Td of the time-dependence has been re- 



lated to the difference Am|i 



m 



I of the 



squared neutrino masses and mi 
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where M m is the mass of the mother ion and 7 = 
1.43 is a Lorentz factor pQ. 

For the calculation of the .EC-decay rate we 
have used the standard weak interaction Hamil- 
ton operator 



(5) 



where U e j are the matrix elements of the mixing 
matrix of massive neutrinos and (x) is defined 
by 



x[^(x)7 M (l- 7 5 )V' e -(2:)] 
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with standard notations [2J. 

The amplitude A(m —> d + v) oi the EC-decay 
m —> d+ v, where m, c? and are the mother ion, 
the daughter ion and a neutrino, has been defined 
as follows 



A(m — > <f + i/) = U e jA(m — > d+ i/j-), 



(7) 



where the coefficients E/ e j testify that the elec- 
tron couples to the electron neutrino. In turn, 
the amplitude A(m — ► d + Vj) is equal to 

A(m d + vj) = - J d 4 x {vjd\H$ '(x)\m) = 

= -(2^6^ (k d + kj - k m ) (vjd\H$(0)\m) = 
= {2^)H {4 \k d + k ] -k m )M (m^d + Vj) (8) 

with M(m -xl+Uj) = -(ujd\H$(0)\m). 

Recently Giunti has criticised this explanation 
[3]. According to Giunti [3], the correct neutrino 
wave function in the final state of the -EC-decays 
of the H-like ions m — > d + v should be taken in 
the form 

\Mt))= ' 



(9) 



(10) 



v 3 

with Ak(t), defined by 

Ak(t) = —i J dr (v k d\Hw(T)\m), 
Jo 

where H\y(t) is the weak interaction Hamilton 
operator Eq.([5]). 

As has been shown in |4J, such a wave func- 
tion contradicts the principles both of standard 
time-dependent perturbation theory [6]-[9] and 
of quantum field theory [1011 lj . 

In the recent paper [5] Giunti has undertaken a 
new attempt to refute the explanation, proposed 
in [2], of the experimental data by GSI pQ. Be- 
low we comment on Giunti's analysis of the EC- 
decay rate. 

The EC decay rate 

According to Giunti [5], the decay probability 
Pm — > d+v) defined by 

Pm^d+v = \A{m^d + v)\ 2 (11) 



with A(m — > d + v) given by Eq.((7j), is not equal 
to 



Pjti — *d-\-v 

^|E/ ej |V(r7wd+^)| 5 



(12) 



In addition Giunti claims that the decay proba- 
bility P m ^,d+u Eq. (fTTj) has an incorrect massless 
limit rrij — > 0, namely 

P m ~,d+u = lira \A(m — > d + v)\ 2 = 



= \A(m^d+u)\ 2 SM ^U ej 

3 

whereas the correct limit is 

P m ~,d+u = lira \A(m — > d + v)Y 

rrij— >0 

= \A(m^d + v)\l M Y,\U ei \ 2 = 



(13) 



\A(m -> d + v) 



SMi 
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where \A(m — > d+ ^OIsmj calculated in the Stan- 
dard Model of electroweak interactions of heavy 
ions, is equal to 

\A(m -» d + v)\%M = 1™ \M m -* d + v i)f- (15) 

The incorrectness of these assertions is clearly 
seen if one takes into account correctly the con- 
tribution of the ^-functions S^(k d + kj — k m ), 
describing the conservation of energy and 3- 
momentum in the EC-decays. 

Substituting Eq.® into Eq.lfTTj) we get 

P m ^d+u = \A(m^d + v)\ 2 = 

= J2 \U e3 \ 2 [(2-K)H {A Hk d + hi - k m )f 

j 

x \M(m -» d + v 3 )\ 2 + 2 Re[U* e U ej 

xM*(m — > d + Vi)M{m — > d+ Vj)] 

x[(2ir) 4 6( 4 \k d + h-k rn )} 

x[(2^)V 4 )(fc d + fc,-fc m )]. (16) 

For subsequent calculations one has to use the 
relations [10] 



[(27r)V 4 )(fc d + fc,-fc m )] 2 = 
= VT(2ir)W>(k d + k j -k m ), 



(17) 
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where VT = {2n) 4 S^ (0) is the volume of the 
space-time [TU], and 



x[(2n) 4 S^(k d + k 3 -k m )} = 



(18) 



The relation Eq. (fT5)) is valid, since fc, ^ fcj. Due 
to the relations Eqs. fTT)) and (JT5]) the decay prob- 
ability Pm^d+v takes the form 

Pm^d+ V = \A(m -^d + u)\ 2 = VTJ2 M 

3 

d 



x (27r) 4 (5 (4) (fc d + kj - k m ) \M(m 
= Y J \Ue ] \ 2 \A(m^d + Vj )\\ 



(19) 



Thus, energy and 3-momentum conservation lead 
to the decay probability, required by Giunti [5]. 



It is obvious that in the limit m,- 



the de- 



cay probability Eq. lfTO]) reduces to the form of 
Eq.fUlD. 



= VT(2ir) 4: S ( - 4 \k d + k v 
x|7W(m^d + l ,)|| M ^|f/ ei | 2 
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where k u = {E v ,k) with E v = \k \ and \M.(m 
Ism 



d + u )\sm is s iven h y 
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The -EC-decay constant is defined by 



X 



EC 
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Thus, the EC-decay constant XeCi determined 
by Eq. (f2"2")) . is equal to the EC-decay constant, 
calculated in [T^] (see also [2]) within the Stan- 
dard Model of the electroweak interactions of 
heavy ions. 



Giunti's wave function of neutrino with 
lepton flavour i 

Now let us make comments on Giunti's wave 
function of neutrino in the final state of the EC- 
decay. In addition to the critique, expounded in 
[4], we would like to emphasize that Giunti's wave 
function of the neutrino vg with the lepton flavour 
t depends on the initial and the final states of the 
reaction Ii — > If + vg in which the neutrino vg is 
produced, where Ii and // are not necessary one- 
particle states. In order to accentuate this point 
we propose to rewrite the wave function Eq.([9]) 
specifying the initial and final states as 



W (t))ii if 



k 



(23) 



with A k (t) ItI given by 



A k (t)=-i / dT (VkIf\H W (T)\Ii, 



(24) 



Hence, the neutrinos (vAli* and {vg)i>j> , pro- 
duced in two different reactions Ii — > If + vi and 
I\ — >■ I'f + i/£, are two different particles. They 
are not stable and the probability of the transi- 
tion {vg)i i i f < — ► {y '1)1^1 ' is equal to 



P{v, 



hi, 



I'J 



V f ) = \i' i i' f (Mt)\Mt))hi f \ = 



A k (t) 



hi, 



hit 



\A f (t) 



(25) 



i'J' 



where Ak(t)i,i f ^ Ak(t)i'i' by definition due to 
different initial and final states of the reactions 
Ii -> If + {v t )hi s and I[ -> I' f + (vg) ri r f , re- 
spectively. Since the number of initial and final 
states {hlf) of the reactions producing neutrinos 
(vg)i i i f with a lepton flavour I is infinite, so, ac- 
cording to Giunti |3l5j . there is an infinite set of 
neutrinos {yg)i i i 1 with a leptonic flavour £. 

Giunti's explanation of "Darmstadt oscil- 
lations" 

According to Giunti [5], the interference term 
in the EC-decay rate m — > d + v e comes from the 
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mixing of the different mass-states of the mother 
ion m. For the wave function of the initial state of 
the mother ion, which is not an eigenstate of the 
mass-operator, Giunti has proposed the following 
expression 

\m) = cos 6\m')+ sm6 \m"), (26) 

where \m!) and \m") are two states of the mother 
ion with masses M m i and M m » , respectively, and 
9 is a mixing angle. This means that the initial 
state of the mother ion is a coherent state of two 
eigenstates of the mass-operator |m') and |m"), 
respectively. The final state of the EC-decay is 
defined by the wave function \d,v e ), where v e is 
a massless electron neutrino. Since Giunti's cal- 
culation has no relation to the real calculation of 
the -EC-decay rate, which is needed for the com- 
parison with the experimental data on the rate of 
the number of daughter ions [2] , below we give a 
calculation of the EC-decay rate within Giunti's 
approach in detail. 

According to standard time-dependent pertur- 
bation theory [6]-[9], the amplitude of the m — * 
d + v e decay is defined by (see also Eq.(fIUj) and 



A(m —> d + v e )(t) = 

= -i dr (is e ,d\H w (T)\m), 
Jo 



(27) 



where \m) is the wave function Eq. (|26|) and the 
weak interaction Hamilton operator Hyy (t) takes 
the form [T2] 

H w (t) = -JLv ud J d 3 x$ n (x)^(l - g A j 5 )Mx)] 



x[i t (^(i-7 5 )4-W]- 



(28) 



Suppose that the wave function of the neutrino 
in the EC-decay is a plane wave. In this case the 
amplitude of the EC-decay is equal to [B"]- |9I12] 



A(m -»• d + i/ e )(t) = - VS y 1 2M m 2E d {q)E VB (k) 

(A§L, 



iAELj 

cosOe 1 2 L 



x — A — ^- +sm6e l ^~ l 



( ae ' 

V 2 



x 6 



(29) 



where AE' = E d (q) + E Ue (k ) - M m , and AE" = 
Ed(q) + E Vc (k ) — M m ii , q and k are 3-momenta 
of the daughter nucleus and the neutrino, E d (q) 
and E Vc {k) are the energies of the daughter ion 
and the neutrino, respectively, A^gt is the nu- 
clear matrix element of the Gamow-Teller transi- 
tion and ) is the wave function of the bound 
electron in the H-like mother ion, averaged over 
the nuclear density [12]. 

The rate of the neutrino spectrum is defined by 

m 

dN Ve (t) If d 3 q 



dt 



2M„ 



2F 



2nf2E d {q) 
l — Y, j t \A{m^d + v e ){tt 



2E + 1 



^e.jmgtI 2 !^)! 2 



I 2 „ sin(AE'i) „ . o „ sin(AE"t) 
x < 2 cos 2 9 —\ + 2 sin 2 6 - ' 



sin 26 



(AE') 



/AE' 



AE' 
2 



V 2 



(AE") 



AE" ) t 



sin 



( AE '' 

V 2 



f AE" 



AE'U 



(30) 



Here we have used the relation 

[(2vr) 3 ^(q + fc)] 2 = V(2ir) 3 5W(q + k), (31) 

where (2n) 3 S^ (0 ) = V is the normalisation vol- 
ume [10]. For sufficiently long time we get [6]-[TTj 

x I cos 2 9 (2ir)6(AE') + sin 2 (2vr) 6{AE") 
+ sin 26 



x cos 



(27r)<5(AE') + (27r)£(AE") 
(AM m t)}, 



(32) 



where AM m = M m , - M m », AE' = E d (k) + 
E Ue (k)-M m , and AE" = E d (k)+E Ve {k)-M m „ . 
The EC-decay rate A^p(i) from the coherent 
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state \m) is defined by 



d 3 fc 1 dJV Ve (t) 



Substituting Eq. (f32f into Eq. (|33f and integrating 
over the neutrino phase volume we get 

X^(t) = A E c{l + 2sin20cos(AM m i)}, (34) 

where A^c is the -EC-decay constant, calculated 
in [12]. 

Thus, one can show that for the initial state 
of the mother ion, given by Eq. ([2^|) . the EC- 
decay rate has a periodic interference term with a 
period Ta, defined by the mass difference AM m = 
M m i — M m n , which is equal to 



AM„ 



T d c 2 



8.45 x l(T ib eV/c 



(35) 



where 7 = 1.43 is the Lorentz factor pQ. A reduc- 
tion of the EC-decay rate Eq. (f34l) to the experi- 
mental shape p] 

/2?rf 



A 



EC 



(t) = \ec{ 



1 + a-EC cos 



■</>)} (36) 



with a ec — 0.20 pQ can be carried out by chang- 
ing \m') -> e^' \mf) and \m") -> e^" \m"), giv- 
ing <j>=tf>' - ((>", and setting 9 ~ 2.87°. 

The problem of such an explanation of the 
"Darmstadt oscillations" is as follows. If there 
exist two states of the mother ion \m!) and 
\m") with a mass-difference AM m — 8.45 x 
10 _16 eV/c 2 , giving the contribution to the EC- 
decay through the coherent state |m), given by 
Eq. (|26|) , the contribution to the EC-decay should 
be also from the coherent state \rh) 

|m) = - sin 1 to') + cos0|to")> (37) 

which is not also an cigenstate of the mass- 
operator and orthogonal to the state \m). The 
coherent state |m), given by Eq. (j37]) . can be 
produced in the system of mother ions on the 
same footing as the coherent state |m), given by 
Eq. (|2"6")) . Indeed, the mother ions, injected into 
the Storage Ring, are produced by means of a fast 
projectile fragmentation with a statistical popula- 
tion of the states \m') and \m"), which are eigen- 
states of the mass-operator. Hence, the probabil- 
ities P m and Pf n of the appearance of the coher- 
ent states I to) and |m), related by P m + Pm = 1, 



should be equal P m — P m = i due to a principle 
indistinguishability of these states. 

The EC-decay rate A^(t) of the EC-decay 
d + v e from the coherent state |to) is equal 



m 
to 

A^ (t) = X E c{ 1 - 2 sin 29 cos(AM m t) | . (38) 

The total EC-decay rate, caused by the EC- 
decays of the H-like heavy ions from the states 
I to) and I to), is defined by 

\ EC (t) = PrA^M*) + P rn ^EC (*) = 

= \ EC {l + 2 sin 29 (P m - P m ) cos(AM ro t)}. (39) 

Since there is no physical reason for P rn 7^ P m , 
setting P m = P^ one gets no interference terms 
in the EC-decay rate of the H-like heavy ion in 
the approach proposed by Giunti [5]. 

Summary 

We have shown that Giunti's critique of our ap- 
proach is based, technically, on the missing of the 
5— functions, responsible for the conservation of 
energy and 3-momentum in the EC-decay and, 
globally, on the misunderstanding of the standard 
procedure for the calculation of the decay rates. 

Giunti's wave functions for neutrinos in the fi- 
nal state of the EC-decay of the H-like heavy 
ions or generally in any weak interaction produc- 
ing or absorbing neutrinos make no sense, since 
they require an infinite number of neutrinos with 
a lepton flavour I. These neutrinos are not stable 
and oscillate with a finite probability. 

As regards Giunti's explanation of the "Darm- 
stadt oscillations" we assert the following. Apart 
from the existence of a superweak interaction, 
leading to the mass-splitting of the H-like heavy 
ions of order O(10~ 15 eV/c 2 ), which is hardly pos- 
sible in reality, in Giunti's approach the total 
EC-decay rate of the H-like ions should be de- 
fined by EC-decays from two coherent states \m) 
and I to), given by Eqs. ([26|) and ([37]) . respectively. 
These coherent states can appear in the system 
of mother ions with probabilities P m and P m , re- 
spectively, and related by P m + P m = 1. The 
equality P m = P m = |, caused by a statistical 
equivalence of the coherent states \m) and \rh) in 
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the system of mother ions, injected into the Stor- 
age Ring, shows the absence of the interference 
term in the EC-decay rate 

\ EC (t) = P m X^ (t) + P m \f£ (t) = X EC - (40) 

Thus, Giunti's explanation of the "Darmstadt os- 
cillations" has no physical ground and it is erro- 
neous by definition. 

Our analysis of Giunti's explanation of the 
"Darmstadt oscillations" can be formulated more 
generally as "a non-existence of periodic time de- 
pendent interference terms in the -EC-decay rates 
of the H-like heavy ions for the coherence in the 
initial state of the mother ions". Indeed, the 
states |m) as well as the orthogonal state |m) can 
be treated as coherent states of \m!) and \m"), 
produced in the statistical system of mother ions 
injected into the Storing Ring. Due to a statis- 
tical equivalence of these states, the probabilities 
P m and P m of the appearance of the coherent 
states \m) and |m) in the system of mother ions 
should be equal P m = Pm. This prohibits the 
appearance of the interference term in the EC- 
decay rate (see Eq.j39]) and Eq.lj3U|l). 

We thank M. Faber, T. Ericson and N. Troit- 
skaya for fruitful discussions. One of us (A. 
Ivanov) is grateful to N. Ivanov for discussions 
of physics of Storage Rings. 
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